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Compulsive behaviors in obsessive-compulsive disorder (OCD) may be related to deficits in reward processing 
mediated by corticostriatal circuitry, a brain network implicated in the pathophysiology of OCD. Performing com- 
pulsive actions can be perceived as a reward to OCD patients because it temporarily reduces the anxiety pro- 
voked by obsessions. Although most OCD literature provides evidence of altered regional activity in these 
corticostriatal circuits, very little is known about the connectivity between individual regions of the 
corticostriatal-limbic circuits, including the cognitive and affective neural circuitry associated with OCD. Thus, 
this study investigated the differences in functional connectivity (FC) patterns in this network during 
resting-state and incentive processing. Nineteen patients with OCD and 18 well-matched healthy controls 
were scanned during resting-state and a monetary incentive delay task (task state). FC was assessed using 
both voxel-wise and region-of-interest (ROI)-wise analyses. Voxel-wise FC analysis with the nucleus accumbens 
seed revealed that patients with OCD exhibited increased FC between the nucleus accumbens and the lateral 
orbitofrontal cortex during resting-state. Additionally, these patients showed decreased FC between the nucleus 
accumbens and limbic areas such as the amygdala during incentive processing. Exploratory ROI-wise FC analysis 
revealed that OCD patients demonstrated enhanced FC between the nucleus accumbens and the lateral 
orbitofrontal cortex and increased total connectivity of the lateral orbitofrontal cortex during resting-state. 
Additionally, patients showed alterations in FC between resting and task state. This study provides evidence 
that patients with OCD have altered FC in the corticostriatal-limbic network, particularly in striatal-amygdala 
and striatal-orbitofrontal circuitry, during incentive processing and resting-state. These findings also emphasize 
that functional connections in the network are modulated by affective/motivational states and further suggest 
that OCD patients may have abnormalities of such modulation in this network. 

© 2013 The Authors. Published by Elsevier Inc. All rights reserved. 



1. Introduction 

One of the striking features of obsessive-compulsive disorder 
(OCD) is the occurrence of repetitive behaviors (compulsions) that 
diminish anxiety-inducing thoughts (obsessions). These compulsions 
produce temporary relief from anxiety that is interpreted by patients 
as a reward, but they never feel fully safe (Cavedini et al„ 2006). 
OCD patients also evidence cognitive and behavioral inflexibility and 
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exhibit deficits in goal-directed and adaptive behaviors, which may 
be seen in alterations of behavioral responses when an outcome or 
environment is altered (Chamberlain et al„ 2007; Gu et al„ 2008; 
Page et al„ 2009); these typical clinical features and maladaptive 
behaviors are commonly thought to stem from an impaired ability to 
process rewards (Cavedini et al., 2006; Huey et al., 2008), because 
reward processing is responsible for appropriate responses to rewards 
or punishments during incentive-based learning and for the develop- 
ment of numerous goal-directed and adaptive behaviors by improving 
future performance and emotion regulation (Haber and Knutson, 
2010). The neural circuits mediating reward processing overlap with 
several brain regions implicated in the pathophysiology of OCD includ- 
ing the basal ganglia (BG), thalamus (THAI), and anterior cingulate cor- 
tex (ACC) as well as areas of the prefrontal cortex (PFC) including the 
orbitofrontal cortex (OFC) and dorsolateral prefrontal cortex (DLPFC) 
(Atmaca et al., 2011; Kwon et al, 2009; Saxena and Rauch, 2000). 
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Particularly, the nucleus accumbens (NA), as part of the ventral striatum 
(VS), constitutes a central region for reward processing (Choi et al., 
2013; Haber and Knutson, 2010) and a primary target of deep brain 
stimulation (DBS) therapy for refractory OCD (Denys et al., 2010). 
Recent meta-analyses of functional imaging using positron emission 
tomography, single photon emission computed tomography (Whiteside 
et al., 2004), and voxel-based morphometry studies (Radua et al., 2010) 
found reliable differences in radiotracer uptake and gray matter volume 
in the VS of OCD patients compared with controls. Abnormalities of 
dopaminergic function, which plays a crucial role in reward processing, 
have also been reported in OCD patients (Denys et al., 2004; van der 
Wee et al., 2004). More recently, abnormal activity during reward pro- 
cessing was identified in the OFC and VS of OCD patients (Chamberlain 
et al., 2008; Figee et al., 2011; Jung et al., 2011; Kaufmann et al., 2013; 
Remijnse et al., 2006). Therefore, an examination of reward processing 
and its neural correlates in OCD patients may be beneficial in understand- 
ing the pathophysiology and treatment of OCD (Huey et al., 2008; Milad 
and Rauch, 2012). 

Converging evidence from functional imaging studies utilizing uni- 
variate general linear model (GLM) analysis reveals altered regional 
activity in corticostriatal circuits and points to their involvement in 
OCD (Gu et al., 2008; Han et al., 2011; Page et al., 2009). However, 
very little is known about the interactions between these regions in 
OCD patients. Recent advances in mapping neural connectivity at the 
circuit level have enabled the exploration of abnormalities in inter- 
regional interaction among corticostriatal regions. These studies gen- 
erally use functional connectivity (FC) analysis that captures the 
temporal structure of resting-state or task-related network dynamics 
(Harrison et al., 2009; Sakai et al., 2011; Stern et al., 2011). Recent 
research has indicated a role for specific brain regions and circuits 
involved in the affective-motivational aspects of behavior of healthy 
controls and OCD patients. These areas include the amygdala (AMY), 
hippocampus (HIP), and insula (Kwon et al., 2009; Menzies et al., 
2008) and anatomically and functionally differential OFC subregions 
including the medial OFC (med-OFC) and the lateral OFC (lat-OFC; 
Milad and Rauch, 2012; O'Doherty et al., 2001 ), as well as the interac- 
tions between these regions, which collectively represent an affec- 
tive and cognitive circuit underlying the pathophysiology of OCD. 
However, as of yet, there has been no specific investigation of func- 
tional interactions between these regions during affective and moti- 
vational processing in OCD. 

A study of FC during a monetary incentive delay (MID) task is 
employed here and may help resolve this issue as the task (i) elicits 
both affective and neural responses to reward and punishment and 
(ii) is associated with the corticostriatal regions implicated in OCD 
including the NA, ACC, anterior insula (A1S), and other limbic regions. 
Previous studies using the MID task in OCD patients reported inconsis- 
tent results in the regional NA activity during reward anticipation- 
decreased activity (Figee et al., 2011) and no group difference 
(Jung et al., 2011 ; Kaufmann et al., 2013). On the other hand, FC anal- 
yses have been shown a consistency in the results of increased ven- 
tral striatal FC in OCD patients during resting-state (Harrison et al., 
2009, 2013; Sakai et al., 2011). Based on these reports, neural dys- 
functions in OCD patients during incentive anticipation may lie in 
functional interactions between brain regions rather than regional 
activity. Therefore, in this study, MID task data (Jung et al., 201 1 ) 
were reanalyzed using FC analysis, specifically a beta-series correla- 
tion technique (Rissman et al., 2004), to examine event-related 
changes in whole-brain voxel-wise FC analysis of the NA. The 
beta-series correlation method can disassociate functional networks 
corresponding to the distinct stages of a cognitive task. 

Recent studies have investigated FC not only during cognitive tasks 
but also during resting-state, and the findings suggested an interac- 
tion between resting-state activity and stimulus-induced activity 
(Northoff et al., 2010). However, few studies have investigated any 
differences in FC between resting-state and task-performance state 



(Arbabshirani et al., 2013). Evaluating the differences in FC between 
these states may clarify how this network is modulated by task 
demands. Therefore, resting-state data acquired from the same sub- 
ject performing the MID task were analyzed to investigate whether 
there are group differences in the behavior-dependent modulation 
of the network. 

The aim of this study is to examine whether OCD patients exhibit 
altered FC in the ventral corticostriatal network, which is relevant to 
both OCD and incentive processing, during incentive task state and 
resting-state as well as any interaction between these states. Based 
on this prevailing model and previous findings from studies of 
OCD, it is hypothesized that, compared with healthy controls, the 
corticostriatal-limbic regions of OCD patients would show increased 
FC during resting-state and reduced FC during incentive processing. 
Specifically, this altered connectivity in OCD patients is expected to 
be found in the OFC and limbic areas, such as the AMY, which have 
been implicated in affective processing. To address this, first, we 
performed whole-brain voxel-wise FC analysis to generate the FC 
map of the NA for both the MID task and resting-state data. We 
also conducted exploratory FC analyses in ROI-wise manner using 
graph theory to examine FC among brain regions that comprise the 
corticostriatal-limbic reward network (Camara et al., 2009; Haber 
and Knutson, 2010; Milad and Rauch, 2012). Then, we compared 
the results from two different analytic methods to show that the 
altered connectivity in OCD patients persists regardless of which 
analytic method is used. 

2. Material and methods 

2. J. Participants 

Twenty patients with OCD were recruited from the outpatient clinic 
at Seoul National University Hospital. Diagnoses and comorbidity were 
established by experienced psychiatrists using the Structured Clinical 
Interview for DSM-IV Axis-I and -II disorders (SCID-I and -II). Twenty 
comparison subjects who were matched for sex, age, and IQ_ with the 
patients were also recruited. All subjects were right-handed. Exclusion 
criteria included lifetime history of psychosis, brain injury, neurological 
illness, alcohol or substance abuse, and learning disabilities. Of all par- 
ticipants, three subjects (two healthy controls, one OCD patient) were 
excluded from FC analysis because of excessive head motion during 
resting-state. The criteria for excessive head motion were translations 
in excess of 1 mm in any directions or rotations in excess of 1 "around 
any axes in 6 head motion parameters. Demographic, clinical, and be- 
havioral data were recorded for the two groups (Table 1 ). At the time 
of the study, 15 OCD patients were drug naive, and four OCD patients 
were drug free for at least 4 weeks. Note that all subjects in this study 
were those included in our previous study using a univariate analysis 
approach (Jung et al., 2011). The institutional review board of Seoul 
National University Hospital approved the study protocol, and all par- 
ticipants signed an informed consent form prior to their participation. 
Participants were paid according to their task performance (see Partic- 
ipants and clinical and behavioral results in the Supplementary 
material). 

2.2. Monetary incentive delay task 

A modified version of the MID task developed by Knutson et al. 
(2001 ) was utilized to invoke anticipation of monetary gains and losses 
in OCD patients and healthy controls (Fig. SI). While performing the 
task, subjects saw one of three cues during each trial, which led them 
to anticipate potential monetary gain, monetary loss, or no conse- 
quence. Participants' monetary outcome depended on their perfor- 
mance during a simple reaction time task at the end of each trial, 
which involved pressing a button while a visual target was briefly 
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Table 1 

Demographic and characteristic features and behavioral data for subjects. 



Controls Patients t p 

(n = 18) (n = 19) 



Demographic data 



Sex (male/female) 


11/7 


12/7 


0.02 b 


0.90 


Age (years) 


24.83 (3.88) 


25.84 (7.15) 


— 0.53 


0.60 


Education (years) 


14.61 (1.30) 


15.16 (3.82) 


— 0.58 


0.57 


IQ score a 


111.11 (10.43) 


108.89 (13.21) 


0.56 


0.58 


BIS 


15.94 (2.62) 


18.63 (3.77) 


— 2.50 


0.02 


BAS total score 


33.06 (7.76) 


30.79 (5.33) 


1.04 


0.31 


BAS reward responsiveness 


13.44 (3.28) 


13.53 (2.37) 


— 0.09 


0.93 


BAS drive 


9.78 (2.62) 


8.89 (2.83) 


0.98 


0.33 


BAS fun seeking 


9.83 (2.38) 


8.37 (2.67) 


1.76 


0.09 


Clinical date 










BDI 


3.UU (4.5b J 


15./9 (11.25J 


— 4.48 


p < 0.01 


BAI 


4.33 (4.69) 


15.63 (14.99) 


-3.06 


p < 0.01 


Y-BOCS total score 




20.63 (6.58) 






Y-BOCS obsession 




10.89 (4.31) 






Y-BOCS compulsion 




9.74 (4.33) 






Age of onset (years) 




18.00 (6.76) 






Duration of illness (years) 




7.82 (5.57) 






Behavioral data 










Total earnings (in won) 


28611.11 (2913.23) 


28631.58 (2832.56) 


-0.02 


0.98 


Reaction time (ms overall) 


230.05 (18.91) 


227.19 (41.27) 


0.27 


0.79 


Hit rate (% overall) 


58.40 (2.86) 


58.25 (1.7) 


0.20 


0.85 


VAS effort for gain (1-10) 


8.26 (1.31) 


8.62 (1.07) 


-0.92 


0.37 


VAS effort to avoid loss (1-10) 


8.82 (0.84) 


8.35 (1.33) 


1.25 


0.22 



Abbreviations: IQ, Intelligence Quotient; BIS, Behavioral Inhibition System; BAS, Behavioral Activation System; BDI, Beck Depression Inventory; BAI, Beck Anxiety Inventory; 
Y-BOCS, Yale-Brown Obsessive Compulsive Scale; VAS, Visual Analog Scale. 
Values are presented as mean (SD). Independent sample t-test was used. 

a Estimated by Korean-Wechsler Adult Intelligence Scale-Revised (K-WAIS-R). 

b x 2 test was used. 



presented (see Supplementary material for a detailed description of the 
trial structure). 

2.3. Data acquisition 

Functional images (T2*-weighted echo-planar imaging, repetition 
time (TR) = 2340 ms, echo time (TE) = 52 ms, flip angle = 90°, 
field of view (FOV) = 22 cm) were acquired using a 1.5-T Siemens 
scanner (Avanto, Germany). Twenty-five contiguous 5.0-mm axial 
slices with 3.44 mm 2 in-plane resolution were obtained. All subjects 
underwent two resting-state scans (i.e., eyes closed (EC) and eyes fix- 
ating on a cross (EF)), followed by the MID tasks. The fMRI scans for 
resting-state and the MID task were performed in a 4.68-min 
(120 volumes) and in two runs of each 9.56 min (490 volumes in 
total) for each subject. In this study, we used only EF data among 
two resting-state scans because it represents more appropriate 
baseline for task-related FC. After functional image scans, we also ac- 
quired structural Tl image for each subject in the same scanner. 
Tl-wighted images were acquired using a 3D magnetization-prepared 
rapid-acquisition gradient echo (MPRAGE) sequence with following 
imaging parameters (TR = 11.6 s, TE = 4.76 ms, FOV = 23 cm, flip 
angle = 15°, 208 slices, 0.45 x 0.45 x 0.9 mm 3 ). 

2.4. Data preprocessing 

The first four time points of resting-state data were discarded to 
avoid the instability of the initial MRI signal. Similarly, the first three 
time points of the MID task data were discarded. Preprocessing steps 
(implemented in SPM8 (http://www.fil.ion.ucl.ac.uk/spm/)), which 
included slice-acquisition timing, motion correction, and spatial nor- 
malization, were conducted according to a standard template provided 
by the Montreal Neurological Institute (MNI). The normalized images 
were smoothed with a 6-mm Gaussian kernel. For the resting-state 
data, linear detrending and band-pass filtering (0.01-0.08 Hz) were 
performed using REST VI .7 (http://www.restfmri.net/; Song et al„ 
2011b). Additionally, head motion parameters, averaged signals from 



white matter and cerebrospinal fluid, and global brain signals were 
regressed out to remove possible spurious variances (Fox et al„ 2005, 
2009). 

2.5. Functional connectivity analysis 

FC analyses were conducted in voxel-wise and ROI-wise manners 
in terms of a priori and exploratory analysis, respectively. Voxel-wise 
FC was calculated using the correlations between NA seeds and the 
beta-series data in all of the other voxels in the whole brain. For 
exploratory ROI-wise FC analyses, the correlations between ROIs 
was conducted on a 1-to-l (edge- wise) and N-to-1 (region- wise) 
basis. The 1-to-l connectivity analysis was defined as the connections 
between each pair of ROIs, and N-to-1 connectivity analysis was 
defined as the sum of each region's connectivity with all other 
regions. A flow chart of data analysis is provided (Fig. 1A). 

2.5.3. Whole-brain voxel-wise connectivity analysis 

FC maps were created for each participant using a beta-series cor- 
relation analysis approach described previously (Rissman et al., 2004) 
for which the validity has been confirmed (Fornito et al„ 2011; 
Rissman et al„ 2004; Ye et al., 2011). For this analysis, a new GLM 
design matrix using separate covariates was constructed to model 
hemodynamic responses in a particular stage (specific cue) for each 
single trial. The parameter estimates (beta values) of cues were 
extracted to form a set of cue-specific beta-series data, which were 
averaged across voxels within the NA seed ROIs defined using Ana- 
tomical Automatic Labeling (AAL) atlas (Tzourio-Mazoyer et al., 
2002), and then correlated with the beta-series of every other voxel 
in the whole brain for each participant using Pearson's correlation. 
The FC map was created separately for the NA seed in each hemi- 
sphere for each task condition. FC maps for resting-state were also 
produced by averaging the time series in the NA seed and then com- 
puting the Pearson's correlation between the averaged seed time 
series and those from each voxel in the whole brain. Therefore, four 
FC maps (i.e., a specific FC map for gain condition, loss condition, 
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Fig. 1. A flowchart of data analysis and region of interest (ROI) used as nodes in the network analysis. (A) Functional connectivity analyses were conducted in both voxel-wise and 
ROl-wise manners for each behavioral state. After task- and resting-state fMRI data preprocessing, we applied the anatomical template for regional parcellation. Voxel-wise functional 
connectivity was estimated by calculating the correlations between the average time series of the nucleus accumbens seed region and the time-series data from all other voxels in 
whole brain. All correlation maps were transformed to the normal distribution by Fisher's z transform for further statistical inference. For ROI-wise functional connectivity computing 
the correlations between 32 ROIs, network functional connectivity was conducted N-to-1 (region-wise) and 1-to-l (edge-wise) basis. From each ROI, we extracted mean beta-series 
to analyze functional connectivity using Person's R correlation coefficient in each behavioral state. N-to-1 connectivity was defined as the sum of each region's connectivity with all 
other regions, while 1-to-l connectivity was defined as the connections between each pair of ROIs. Finally, between-group comparisons were performed. (B) Schematic representation 
of the 32 regions comprising the corticostriatal reward network considered in the functional connectivity analysis. Thirty-two redefined regions are shown on (C) a semi-transparent 3D 
rendering of an MRI scan of the brain and on ( D) a single subject Tl data. 



neutral condition, and resting-state) were made for the NA seed in 
each hemisphere for each subject. All correlation maps were trans- 
formed to the normal distribution by Fisher's z transform. To make 
FC maps for the gain (loss) anticipation contrast, FC maps of the gain 
(loss) condition were subtracted from those of the neutral condition. 



Group-level analyses were carried out using a random-effects model 
implemented in SPM8. A 2 x 2 full-factorial design [group (controls 
versus patients) by hemisphere (right and left)] was performed for 
each behavioral state separately. Between-group comparisons were re- 
stricted to (i.e., were masked by) voxels showing significant connectivity 
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with the NA seed within either controls or patients depending on 
the specific behavioral state as a mask image. The mask image was 
thresholded at p < 0.001, and a topological cluster-size false-discovery 
rate (FDR) of p < 0.05 was then used to control for multiple comparisons 
across the whole brain (Chumbley and Friston, 2009; Chumbley et al., 
2010; Fig. 2). Significant areas for between-group effects were reported 
by thresholding individual voxels at p < 0.005 (uncorrected) and then 
applying a subsequent cluster-size threshold based on Monte Carlo sim- 
ulations using AlphaSim, which resulted in a corrected threshold of 
p < 0.01 (Fig. 3). To investigate changes in FC of the NA between 
resting-state and task-performance state in the two groups, a 2 x 2 
full-factorial design was also used for the FC difference map obtained 
by subtracting the gain (loss) anticipation map from the resting-state 
map for two groups. This analysis was limited to only those regions 
showing significant connectivity with the NA within either controls or 
patients during resting-state. 



2.5.2. Exploratory ROI-wise connectivity analysis 

2.5.2.1. Regions of interest. Based on literature investigating the func- 
tional neuroanatomy of incentive processing (Camara et al., 2009; 
Haber and Knutson, 2010; Milad and Rauch, 2012), the following 
brain regions belonging to the corticostriatal-limbic reward network 
were considered to be of particular interest: NA, AMY, HIP, substantia 
nigra (SN), AIS, caudate nucleus (CAU), putamen (PUT), globus pallidus 
(CP), Brodmann's area (BA) 9 (BA9) of the DLPFC, frontopolar cortex 
(FP, BA10), med-OFC (BA11), ventral ACC (BA24), dorsal ACC (BA32), 
BA 46 of DLPFC, lat-OFC (BA47), and thalamus (THAL). Thus, 16 seed 
regions were defined in each hemisphere (32 seed regions in total). 
To improve functional specificity, these ROls were defined based on 
the AAL atlas for subcortical areas (Tzourio-Mazoyer et al., 2002) and 
the BA atlas for the frontal areas, except for the AIS which was manually 
refined from a previous study (Song et al., 2011a; Fig. 1B-D). 




Fig. 2. Significant within-group functional connectivity maps of the nucleus accumbens seed region during (A) resting-state, (B) gain anticipation, and (C) loss anticipation. Red 
indicates positive correlation maps in healthy controls; green, positive correlation maps in patients with OCD; yellow, the spatial overlap of positive correlation maps between 
two groups; blue, negative correlation maps in healthy controls; cyan, negative correlation maps in patients with OCD; L, left hemisphere; R, right hemisphere. Results are displayed 
at p < 0.05, topological FDR corrected. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.) 
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Fig. 3. Significant between-group differences in ventral striatum functional connectivity. This figure indicates significant group differences in voxel-wise functional connectivity of the 
nucleus accumbens seed during (A) resting-state. (B) gain anticipation and (C) loss anticipation. Results are displayed at a corrected p < 0.01, determined by Monte Carlo simulation 
using AlphaSim. 



2.5.2.2. 1-ro-I connectivity analysis. After extracting the 32 ROIs for 
each subject, the FC between each pair of the 32 ROIs was computed 
using Pearson's correlation. A 32 x 32 matrix was obtained for each 
subject and each condition, with each element representing the 
strength of the FC between the two corresponding brain regions within 
the network. The resulting correlation coefficients were then normal- 
ized to z-scores with Fisher's r-to-z transformation. A network-based 
statistic (NBS) approach, as proposed by Zalesky et al. (2010a), was 
used to identify impaired sub-networks in patients for each behavioral 
state and FC differences between behavioral states (see Supplementary 
material). 

2.5.2.3. N-to-1 connectivity analysis. The analysis method described by 
Jiang et al. (2004) was applied. Based on graph theory, all the ROIs 
were referred to as "nodes" and the connections among them were con- 
sidered the "links" in the network. The total connectivity degree T, of a 
node i in a graph was the sum of all the connectivity degrees between 
i and all other nodes. This describes the amount of information received 
by node i from a particular network, specifically, the corticostriatal- 
limbic network in this study. The T value of each ROl was calculated 
and further normalized by dividing the sum of degrees of all considered 
nodes in the network, denoted by T. The differences inF for each ROI be- 
tween the two groups for each behavioral state were tested using 
two-sample t-tests. Additionally, two-sample r-tests were used to com- 
pare the T difference (task minus rest) in the two groups. The signifi- 
cance threshold for between-group differences was set at p < 0.05 
and corrected for multiple comparisons based on permutation testing 
using Ptest software with 10,000 permutations (Camargo et al., 2008). 
The procedures and corresponding mathematical explanations for the 
analysis are provided in the Supplementary material. 

3. Results 



and BG ( Fig. 2A). The spatial distribution of these brain regions is consis- 
tent with the so-called reward system (Haber and Knutson, 2010). 
Compared with patients, healthy controls demonstrated greater FC of 
the NA to the right PUT, parahippocampus, and left posterior insula, 
whereas OCD patients demonstrated greater FC of the NA to the right 
lat-OFC (BA47; Fig. 3A). A number of significant group differences for 
each behavioral state were found in several brain regions (Table 2). 
During gain anticipation in both groups, the NA showed significant FC 
with the AMY, HIP, THAL, insula, BG, and cerebellum (Fig. 2B). Addition- 
ally, OCD patients showed a pattern of significantly increased FC of the 
NA with the posterior insula and occipital regions. Healthy controls had 
greater FC of the NA with the left AMY positioned adjacent to the AIS, 
middle frontal cortex, and midbrain. OCD patients had greater FC of 
the NA to the posterior insula, occipital cortex, and the cerebellum 
(Fig. 3B). During loss anticipation, the NA of healthy controls was func- 
tionally connected to brain regions similar to the FC map during gain 
anticipation, whereas OCD patients showed a set of regions that mostly 
did not overlap those of healthy controls, including the striatum, THAL, 
and cerebellum (Fig. 2C). Healthy controls had greater FC of the NA to 
the bilateral AMY positioned adjacent to the AIS, left mediodorsal 
THAL, middle temporal cortex, and the cerebellum, whereas OCD pa- 
tients had greater FC of the NA to the occipital cortex (Fig. 3C). When 
comparing two groups on FC difference maps between resting-state 
and loss anticipation, significant group differences were found in FC of 
the NA to the med-OFC and rostral ACC (Fig. 4). 

Networks that were negatively correlated with the NA were also 
assessed for both groups. During resting-state, only healthy controls 
exhibited a negative correlation between the NA and regions consti- 
tuting the task-positive network (TPN; Fig. 2A; Fig. S2). Compared 
with controls, OCD patients showed increased FC in TPN regions 
(Fig. S2). There were no significant negative correlations with the 
NA for other conditions in either group. 



3.1. Voxel-wise connectivity 



3.2. Brain-behavior associations 



During resting-state in both groups, the NA was functionally 
connected in a positive manner to several brain regions including the 
bilateral OFC, AIS, AMY, HIP, THAL, SN, ACC, middle cingulate cortex, 



Voxel-wise linear regression analyses in SPM8 were performed to 
assess the extent to which patients' overall symptom severity (total 
Y-BOCS score) may be associated with the FC strength of the NA for 
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Table 2 

Regions demonstrating significant differences between the two groups in voxel-wise 
functional connectivity with the nucleus accumbens seed in each behavioral state. 

Region L/R MNI coordinate 

x.y.z t z 



Resting-state 
Positive network 
Controls > patients 



Putamen 


R 


27, —21, 3 


5.05 


4.6 


Parahippocampus 


R 


21, -18, -18 


3.77 


3.59 


Posterior insula 


L 


— 39, 0,-6 


3.58 


3.42 


Controls < patients 










Cerebellum 


L 


-21, -81, -54 


4.58 


4.27 


Orbitofrontal cortex 


R 


27, 30, -9 


4.29 


4.03 


Gain anticipation 










Controls > patients 










Amygdala 


L 


-33, 9, -15 


4.27 


4.01 


Middle frontal cortex 


L 


-27, 6, 45 


3.63 


3.46 


Midbrain (Substantia nigra) 


L 


0, —18, -12 


3.29 


3.16 


Controls < patients 










Calcarine 


R 


24, -69, 15 


5.07 


4.66 


Cerebellum 


R 


15, —66, 3 


4.48 


4.18 


Calcarine 


L 


— 6, — 54, 3 


4.10 


3.87 


Posterior insula 


L 


-39, 0, 12 


4.06 


3.83 


Cerebellum 


L 


-15, -51, -12 


3.74 


3.56 


Posterior insula 


R 


45, -6, 6 


3.53 


3.38 


Temporal pole 


L 


-39, 0, -18 


3.35 


3.22 


Loss anticipation 










Controls > patients 










Globus pallidus. Amygdala 


R 


30, -18, -6 


4.70 


4.37 


Middle temporal cortex 


R 


42, -6,-21 


4.55 


4.24 


Thalamus 


L 


-3, -24,0 


4.53 


4.23 


Middle temporal cortex 


L 


-42, -30, 0 


4.33 


4.06 


Thalamus 


L 


-3, 0, -3 


4.21 


3.96 


Cerebellum 


R 


48, -57, -36 


3.85 


3.66 


Inferior frontal cortex 


L 


-48, 18, 8 


3.78 


3.59 


Fusiform 


R 


42, -27, -21 


3.68 


3.51 


Cerebellum 


L 


-42, -66, -27 


3.67 


3.50 


Amygdala 


L 


-33, 6, -15 


3.67 


3.50 


Controls < patients 










Calcarine 


R 


24, -69, 15 


4.09 


3.86 



Abbreviations: L, left; R, right. 



each state. Correlations were considered significant if they exceeded 
p < 0.001 uncorrected and reflected a cluster size of at least 18 voxels 
atp < 0.01 (AlphaSim corrected). During resting-state, OCD patients' 
overall symptom severity was positively correlated with FC between 
the NA and med-OFC (x, y, z = 0, 48, -27; z-score, 3.80) and nega- 
tively correlated with FC between the NA and lat-OFC (x, y, z = 21, 
21,-15; z-score, 5.39; Fig. 5). 

3.3. Exploratory ROI-wise connectivity 

Analysis using 1-to-l connectivity in the reward network model 
found significant group differences in FC between several region 
pairs (Fig. 6A-C; Table SI ). While OCD patients exhibited significantly 
increased connectivity in several region pairs relative to controls, par- 
ticularly in FC between PFC regions and BG during resting-state and 
between the striatum and limbic areas during gain anticipation, they 
showed decreased connectivity between the PFC and caudate during 
loss anticipation (Fig. 6A-C). The differences inF among ROls between 
the two groups were provided using N-to-1 connectivity analysis 
(Fig. 6D; Table S2). While OCD patients had significant increases in T 
value of the bilateral lat-OFC during resting-state and of the 
right SN during gain anticipation, they showed significant declines 
in T value of the left NA and right med-OFC (BA1 1 ) during gain antic- 
ipation and of the right DLPFC (BA9) during loss anticipation. Correla- 
tion analyses were also conducted on connectivity values showing 



significant group differences among all analyses and patients' overall 
symptom severity using SPSS. However, there were no significant asso- 
ciations between connectivity measures and symptom severity. More de- 
tailed results and discussions for exploratory ROI-wise connectivity 
analysis were provided in the Supplementary material. 

4. Discussion 

This is the first study to investigate changes in the FC of the 
corticostriatal-limbic reward network, a system relevant to both the 
pathophysiology of OCD and incentive processing in unmedicated pa- 
tients with OCD, during the MID task and resting-state. The main 
findings of this study were as follows: (i) OCD patients had signifi- 
cantly increased FC of NA with the lat-OFC and significantly decreased 
FC of NA with the insula during resting-state; (ii) patients' overall 
symptom severity correlated with the strength of FC between the 
NA and OFC regions during resting-state; (iii) during incentive pro- 
cessing, OCD patients showed decreased FC of NA with the AMY; 
(iv) OCD patients showed aberrant NA connectivity modulation 
between rest and task, particularly in FC of the NA with the med- 
OFC and rostral ACC; (v) consistent with the findings from voxel- 
wise FC analysis, exploratory ROI-wise analyses also revealed signifi- 
cantly increased FC between the NA and lat-OFC in patients during 
resting-state. 

4.1. Ventral striatal functional connectivity 

During resting-state, the pattern of voxel-wise FC to the NA seed 
was quite similar to the reward network reported in previous studies 
(Camara et al., 2009; Haber and Knutson, 2010; Milad and Rauch, 
2012). The current findings reveal that OCD patients had increased 
FC between the NA and lat-OFC compared with controls. A prepon- 
derance of evidence for structural and functional abnormalities in 
this region in OCD patients (Kwon et al, 2009; Saxena and Rauch, 

2000) suggests that the OFC is integral to the neurobiological model 
of this disorder (Chamberlain et al., 2008; Saxena et al., 2001). Addi- 
tionally, this region is critical for adaptive behavior. However, previous 
studies have generally associated the OFC with OCD pathophysiology 
without an explicit distinction between the differential roles of the 
med-OFC and lat-OFC. This is an important distinction, as recent neuro- 
imaging studies have demonstrated functional subdivisions of the OFC 
(Kahnt et al., 2012; Kringelbach and Rolls, 2004). The lat-OFC is thought 
to process negative valence (punishment) and behavioral responses, 
particularly behavioral inhibition in a motivational context (Ursu et 
al., 2008), whereas the med-OFC is thought to process positive valence 
(reward) and emotional processing (Elliott et al., 2000; O'Doherty et al., 

2001) . 

The current findings show that during resting-state, the FC 
strength of the NA to the med-OFC was positively correlated with 
patients' symptom severity, whereas the FC of the NA to a more later- 
ally placed area of the OFC was inversely associated with symptom 
severity. Based on the distinct roles of the two OFC subregions, the 
FC strength between the NA and lat-OFC in OCD patients during 
resting-state may be associated with processing related to the inhibi- 
tion of their clinical symptoms and/or negative affective states associ- 
ated with these symptoms. FC strength between the NA and med-OFC 
may be related to affective states associated with the clinical symp- 
toms. Previous studies showing that larger changes in OCD scores 
are correlated with lower OFC activity during symptom provocation 
(Adler et al., 2000; Rauch et al., 1994) suggest that the OFC may 
play an inhibitory role in OCD. Using resting-state FC, Harrison et al. 
(2009) recently reported greater FC of the NA with both the anter- 
olateral (BA47) and med-OFC in OCD patients. Particularly, dysfunction 
in the lat-OFC is considered a vulnerability marker of OCD, as such activ- 
ity has been identified in unaffected relatives (Chamberlain et al., 2008). 
Therefore, although it is suggested here that the lat-OFC plays a crucial 
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Fig. 4. Regions showing significant group differences on FC difference maps between rest and task. When comparing two groups on FC difference maps between rest and loss 
anticipation, significant group differences were found in functional connectivity of the NA to the med-OFC (x, y, z = —15. 39, -9; t-/z-value = 4.09/3.54) and rostral 
ACC (x, y, z = — 9, 39, 6; t-/z-value = 3.54/3.38). Panel B indicates the connectivity strength of the NA to the med-OFC and rostral ACC during loss anticipation and resting-state in 
each group. We extracted the connectivity strength of the NA to these two regions for each hemispheric seed ROI and each condition, and then these values were averaged for left 
and right hemispheres for each subject. OCD patients showed less modulation of FC between these two states than controls. Error bars represent standard error from the mean. 



role in the pathophysiology of OCD, further studies should clarify 
the specific individual roles of the lat- and med-OFC in subgroups of 
patients with high and low levels of OCD symptoms. 



Recently, researchers have emphasized the central role of brain 
regions such as the AMY, HIP, and insula, which are associated with 
affective/motivational processing in corticostriatal circuits, in the 



A 




Lateral orbitofrontal cortex 




Medial orbitofrontal cortex 
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Total Y-BOCS score 




10 15 20 25 30 
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Fig. 5. Correlations between the functional connectivity strength of the nucleus accumbens to the orbitofrontal cortex and patients' total Y-BOCS scores (p < 0.001) during 
resting-state. (A) During resting-state, patient's overall symptom severity was correlated positively with the FC between the NA and med-OFC (warm color) and negatively 
with the FC between the NA and lat-OFC (cold color) respectively. (B) We extracted the connectivity strength between the nucleus accumbens and OFC regions showing such 
correlations using Marsbar and then assessed Pearson's correlation between the strength and Y-BOCS scores using SPSS. (For interpretation of the references to color in this figure 
legend, the reader is referred to the web version of this article.) 
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Fig. 6. Results from exploratory ROl-wise connectivity analyses. Top figures indicate significant group differences in region-wise (1-to-l) functional connectivity during (A) resting-state, 
(B) gain anticipation, and (C) loss anticipation. To localize specific pairs of brain regions in which functional connectivity was altered in patients, we used a network-based-statistic 
approach (NBS). Bottom panel (D) indicates significant group differences in total (N-to-1 ) functional connectivity degree of each brain region for each behavioral state. To determine 
significant differences in the two groups, we performed permutation-based multiple comparison correction (p < 0.05). R, right; L, left; med, medial; lat, lateral; BA, Brodmann area; 
NA, nucleus accumbens; SN, substantia nigra; OFC, orbitofrontal cortex; DLPFC, dorsolateral prefrontal cortex. AMY, amygdala; HIP, hippocampus; AIS, anterior insula; CAU, caudate 
nucleus; PUT, putamen; GP, globus pallidus; THAL, thalamus. Please see Fig. 1 for abbreviations. 



pathophysiology of OCD (Kwon et al., 2009; Menzies et al., 2008; 
Milad and Rauch, 2012). During gain and loss anticipation states, the 
NA in healthy controls was functionally connected to large parts of 
the affective network (i.e., limbic system) including the AMY, HIP, 
SN, THAL, and BG, whereas FC maps of OCD patients were quite differ- 
ent from those of controls. OCD patients showed increased FC between 
the NA and occipital regions but decreased FC between the NA and 
AMY during both gain and loss anticipation. Neural responses in 
sensory regions, including the occipital cortex, are modulated by 
attention (O'Connor et al., 2002), and it is thought that OCD patients 
focus their attention on incentive cues. An abnormal amygdala re- 
sponse to emotional stimuli was reported in OCD patients, albeit 
with inconsistent results, showing both hyper- (Cardoner et al., 
2011) and hypo-responsiveness (Britton et al., 2010). Consistent 
with a previous study showing less amygdala activation in OCD pa- 
tients in response to multiple emotional stimuli (Britton et al., 
2010), subjects in the present study exhibited decreased FC in the 
AMY to both incentive cues. OCD patients also showed increased FC 
between the NA and posterior insula during gain anticipation, but 
the FC between these regions was decreased during resting-state. 
The posterior insula is known to be a primary interoceptive brain 
region, and it plays an important role in the sense of the physiological 
condition of the entire body (Craig, 2002). In this regard, aberrant FC 
in the posterior insula may be involved in deficits of homeostatic 
processes according to behavioral states. 

A negatively correlated (anti-correlation) network with the NA seed 
was observed in healthy controls only during resting-state (Fig. 2); this 
closely matches the TPN, which is anti-correlated with the DMN. 



However, the anti-correlation network of the OCD group, including 
the occipital region, cerebellum, and THAL, generally did not overlap 
that of controls (Fig. 2; Fig. S2). These differences in anti-correlation 
network patterns resulted in significant group differences in TPN re- 
gions when comparing the FC difference maps between task and rest 
in the two groups (Fig. S2). However, several researchers have recom- 
mended caution during interpretation of anti-correlations when using 
global signal regression (Chang and Glover, 2009; Murphy et al., 
2009) because there has been the debate about enhanced anti- 
correlation networks following global signal regression (Fox et al., 
2009; Murphy et al., 2009). Therefore, further studies should confirm 
this finding in the anti-correlation network using other noise- 
reduction methods such as the component base noise-reduction 
method, which does not rely on global signal regression (Chai et al., 
2012). 

4.2. Methodological considerations in ventral striatal functional 
connectivity studies 

There have been several studies demonstrating changes of ventral 
striatal FC during resting-state in OCD (Fitzgerald et al., 2011; 
Harrison et al., 2009, 2013; Sakai et al., 2011). Harrison et al. (2009) 
reported heightened connectivity of the inferior VS (VSi), corre- 
sponding to the NA, with both the anterolateral (BA 47) and medial 
portion of OFC in medicated OCD patients, while Sakai et al. (2011) 
demonstrated increased FC of the superior VS (VSs), corresponding 
to the ventral caudate, with the med-OFC region in non-medicated 
OCD patients. Fitzgerald et al. (2011) failed to show increased FC of 
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the NA with the OFC in patients. Discrepancies and similarities 
between previous findings and our results may stem from differences 
in resting-state condition, ROI definition, and characteristics of the 
sample (i.e., symptom severity and symptom dimension). First, the 
strength of FC can be affected by resting-state conditions (i.e., eyes 
open (EO), EC, or EF) (Patriat et al, 2013; Yan et al., 2009). While 
three of the previous studies used EC resting-state data (Harrison et 
al., 2009, 2013; Sakai et al., 2011), Fitzgerald et al. (2011) and this 
study used EF resting-state data. Second, the spatial differences of 
seeds (e.g., seed size or location) can also influence the FC strength 
(Cole et al., 2010; Wang et al., 2009; Zalesky et al., 2010b). While 
the NA seed ROIs in this study was anatomically defined using 
AAL atlas, other studies used spherical seed ROIs as described by 
Di Martino et al. (2008). Third, differences in symptom severity and 
symptom dimensions in patients between the present and previous 
studies may have resulted in different strengths of FC. A significantly 
greater FC between the VS and anterolateral OFC regions was found in 
patients with moderate symptoms (the mean total Y-BOCS score < 23, 
Harrison et al., 2009, 2013), while it wasn't observed in patients with 
severe symptoms (the mean total Y-BOCS score > 24, Sakai et al., 
2011). Interestingly, this idea was supported by our findings from 
the correlation analysis; FC strength of the NA to the med-OFC was 
positively correlated with total Y-BOCS score, whereas FC strength of 
the NA to the lat-OFC was negatively correlated with total Y-BOCS 
score during resting-state. Lastly, small sample size, heterogeneity of 
symptom dimensions, and effects of medication in patients may 
have caused limited power to detect the difference in the FC maps of 
the NA between patients and controls in previous studies. To address 
this issue, Harrison et al. (2013) recently investigated the influence 
of major OCD symptom dimensions on FC maps of VSi and dorsal cau- 
date seeds in a large sample of medicated OCD patients (n = 74) and 
controls (n = 74). Similar to their past findings (Harrison etal., 2009), 
they found increased FC of the NA to the anterolateral OFC and 
med-OFC regions in OCD. They also reported reduced FC of the NA to 
the insula in OCD. Our findings are consistent with these results, 
showing increased FC of the NA to the anterolateral OFC and decreased 
FC of the NA to the insula during resting-state. Harrison et al. (2013) 
also demonstrated distinct neuroanatomical relationships between 
certain symptom dimensions and the FC of corticostriatal systems. 
Taken with the aforementioned results and our current findings, the 
increased FC of the VS to the OFC regions may be a common connectiv- 
ity alteration in OCD regardless of the effects of particular symptom 
dimensions (Harrison et al., 2013) and medication. 



4.3. Modulation of functional connectivity during resting and task states 

We investigated how the resting-state FC of the NA changed in 
response to task demands in the two groups. Significant group differ- 
ences on FC difference maps between rest and task were found in FC 
of the NA to the med-OFC and rostral ACC on FC difference maps 
between resting-state and loss anticipation. When extracting the 
connectivity strength within these regions for each condition 
(Fig. 4B), OCD patients, unlike controls, had less modulation of FC 
between these regions. This suggests that functional interactions 
of these regions persist between resting-state and loss anticipation 
(task-performance state). These two regions are involved in emo- 
tional processing and are part of the DMN. Some studies using 
resting-state FC reported positive functional connections between 
the NA and regions of the DMN (Duan et al., 2012; Kelly et al., 
2009). Previous functional imaging studies using GLM analysis 
reported that OCD patients showed reduced DMN deactivation 
during task performance, particularly in the med-OFC (Fitzgerald 
et al., 2010; Stern et al., 2011). Thus, our current finding may be re- 
lated to results showing FC deficits in the DMN of OCD patients 
(Jang et al., 2010). 



4.4. Connectivity between regions belonging to the corticostriatal-limbic 
network as exploratory ROl-wise analysis 

In 1-to-l connectivity analysis, OCD patients show increased FC in 
regional pairs between the BG and PFC regions during resting-state 
and between the striatum and limbic regions during gain anticipa- 
tion. Particularly, increased FC between the NA and lat-OFC (BA47) 
during resting-state was found in OCD patients, which is consistent 
with results from voxel-wise FC analyses. The between-group com- 
parisons in N-to-1 connectivity showed an increased T value for the 
lat-OFC during resting-state and declines in T values in the NA and 
med-OFC, main areas associated with reward processing, during 
gain anticipation in OCD patients. Interestingly, abnormalities in FC 
between the NA and lat-OFC in patients were found in both voxel- 
wise and exploratory ROI-wise FC analyses. Taken together, this 
result suggests that altered FC of the lat-OFC with other regions with- 
in corticostriatal-limbic network is a reliable pathophysiology of OCD 
that is consistently found regardless of which analytic method is used. 

4.5. Limitations 

There are some limitations in this study. First, although patients 
who had comorbid depressive disorder were excluded from this 
study, depression symptoms in patients could have confounding 
effects. However, our main findings remained unchanged when the 
data were reanalyzed with the inclusion of BDI as a covariate, show- 
ing altered FC of the NA with the lat-OFC and insula during resting- 
state and with the AMY during incentive processing. Second, small 
sample size and heterogeneity of symptom dimensions could have 
influenced the present findings. However, our findings were consis- 
tent with results from previous studies that examined the effects 
of symptom dimensions on the FC in a large OCD sample, showing 
increased FC between the NA and OFC regions and reduced FC be- 
tween the NA and insula during resting-state (Harrison et al., 2013). 
Finally, although using seed ROIs defined by brain atlases is prevalent, 
spatial differences between seed ROIs (i.e., seed size) for exploratory 
ROI-wise connectivity analysis could have impacts on FC (Cole et al., 
2010; Wang et al., 2009; Zalesky et al., 2010b). 

5. Conclusions 

Our results reveal specific abnormalities in FC of the NA in OCD, 
showing increased FC of the NA to the lat-OFC during resting-state 
and decreased FC of the NA to the AMY during incentive processing. 
These findings provide direct evidence of altered FC in corticostriatal- 
limbic network during incentive processing and resting-state, particu- 
larly deficits in striatal-amygdala and striatal-orbitofrontal interactions 
respectively. Further, our findings suggest that OCD patients may have 
deficient modulation of FC within the corticostriatal-limbic network 
associated with affective and motivational states. These connectivity 
abnormalities may underpin some of affective and motivational 
disturbances observed in OCD. 
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